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ABSTRACT 


Crack growth inspection data gathered during the service life of the 
C-130 "Hercules" airplane is used in conjunction with a crack propagation rule 
to estimate the distribution of crack initiation times and of initial crack 
sizes. 

Because at early inspections many small cracks are missed a Bayesian 
statistical approach is used to calculate the fraction of undetected initia- 
tion times as a function of the inspection time and the reliability of the 
inspection procedure used. 
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Section I 


INTRODUCTION 

The Lockheed C-130 aircraft has been in use by the U.S. Air Force 
for more than two decades in various theaters of operation under 
several commands. 

The center wing box of the airplane has been periodically inspected 
and the size of cracks growing around rivet holes in the skin has been 
recorded. Small cracks were permitted to grow through several inspec- 
tions. Eventually cracks were repaired or in some cases the whole 
center wing box was replaced. 

The aim of this paper is to reduce the C-130 service-flight inspec- 
tion data by means of fracture mechanics and statistical analysis into 
a form from which the probability of time to crack initiation and the 
distribution of initial flaw sizes for various locations on the aircraft 
can be determined. 

The importance of the time to crack initiation for the purposes of 
reliability analysis and maintenance schedules for aircraft has been 
indicated by numerous authors [1-8]. Though such data has been avail- 
able for laboratory test specimens [9-10] they have not been computed 
from actual service inspection records for a large number of full 
scale structures. For the purposes of crack propagation analysis the 
time to crack initiation, defined here as the time required for a 
crack to grow to an inspectable size, is not needed if the initial 
crack size is known [11-15]. 

Because initial crack sizes are noninspeetable they are computed 
from crack propagation relations by backward extrapolation [11,12]. 
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Previous attempts to define the distribution of time to crack ini- 
tiation for the C-130 and C-141 aircraft [16-18] have not yielded satis- 
factory results. Impel! izzeri [19] discussed methods to evaluate time 
to crack initiation for the F-4 aircraft but the lack of sufficient 
data made the problem difficult. In these reports the data was not 
normalized, that is, the time at first crack sighting was used as the 
time of crack initiation regardless of crack size. This type of 
approach is reasonable only if inspection intervals are short so that 
all the cracks are small and of approximately the same size at first 
sighting. It is, however, impractical and uneconomical to perform in- 
spections at such short service intervals [5,6] . 

The inspection records do not contain cracks smaller than .03 
inches (.0762 cm) in length. Therefore, this size has been arbitrarily 
chosen as the minimum inspection size (crack initiation size) (a^ ). 

Since the objectives of this investigation are to determine the 
distributions of the times of growth to minimum inspection size (crack 
initiation) and of initial flaw sizes, crack sizes are normalized with 
the use of fracture mechanics. A growing crack is extrapolated backwards 
from its first recorded crack size to .03 inches (minimum inspection 
size) yielding the time to crack initiation (t. ) [6]. In a similar 
fashion, the crack is grown backwards to a time equal to zero to yield an 
initial flaw size (a Q ). 

All of the C-130 crack inspection data examined refer to rivet holes 
in the center wing box section of the aircraft (Fig. 1). The center 
wing box was divided into nearly one hundred inspection locations which 
were symmetric about the center line of the aircraft. 
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The original inspection data was received from Warner Robbins Air 
Logistic Center. (ALC/ACDCJ), on a nine track magnetic tape. 

The following information is available for each aircraft: 

aircraft series 

aircraft serial number 

total flying hours (at time of inspection) 

date of inspection 

military command 

military base 

facility where inspected 

number of inspection 

crack locations 

crack sizes 

crack numbers (there may be more than one crack at a 
location) 

In order to perform the discussed computations , crack growth con- 
stants were computed for each individual location and were analyzed 
statistically. 

Two approaches were used to describe the rate of crack growth at 
various locations. One is the more convnon practice of fixing the crack 
growth parameters. The second approach is to let the crack growth 
parameters vary from location to location even though the material, 
geometry and type of loading are similar between locations. 

Only cracks that grew were utilized in the present analysis. These 
comprise approximately fifteen percent of the total number of rivet 
holes and consequently conclusions reached are valid for this reduced 
population. 

Statistical distributions were fitted to time to crack initiation 
and to initial flaw sizes in order to describe their expected behavior 
at each location. Since only the cracked population was considered, 
probabilities are conditional on the presence of a crack [20], 
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The influence of the modes of operation within the various commands 
is reflected in the rates of crack propagation indicating more or less 
severe usage. The data has, as a consequence, been also used to esti- 
mate the parameters of the load distributions applied to the airplanes. 
Because the method of inspection is not a hundred percent reliable 
some cracks present in the structure were presumably not discovered. 
Utilizing an assumed inspection reliability curve combined with a 
Bayes' approach, the distribution of crack initiation times undetected 
during an inspection has been estimated [21,22], 
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Section II 


METHODS OF ANALYSIS 

1 . Time to Crack Initiation and Initial Flaw Size 

Inspection of the aircraft took place at irregular intervals. As a 
result raw inspection data points could not be directly compared because 
the length of time during which cracks have propagated differed for each 
crack. 

Once the crack propagation parameters have been calculated, cracks 
could be mathematically grown to either a common size or to a common 
inspection time. Such a "normalization" procedure permits the determina- 
tion of statistical distributions for "crack initiation times" defined 
here as the time at which the crack size was .03 in (minimum inspectable 
crack length), and for "initial crack sizes", the size of cracks at time 
equal to zero (Fig. 2). 

In order to determine the time to crack initiation (t.) and the 
initial flaw size (a Q ) the "Power Law" crack growth relation [23] 

da/dN = C(Ak) n (1) 

was used. 

Here Ak = AaVira f(l/r) = Stress Intensity Range (2) 

da/dN = crack growth per cycle 
C & n = material constants 
A a = stress range 
a = crack length 

f O/r) = geometric correction factor 
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Eq. 1 may be rewritten utlizing Eq. 2 and a rate-of-cycl ing term, 
dN/dt as 


/ \ n 



da/dt = C ( Aa/ira" f(l/r) J 

dN/dt 

(3) 

/ Y n 



= C ( Aa/iT f(l/r)j 

dN/dt a n/2 


/ \ n 



Letting C* = C f Aa/ir f(l/r)J dN/dt 


(4) 

and 



n* - n/2 


(5) 

where C* is in units of [in./in. n ^ 2 /fl ight hours], Eq. 3 becomes 


da/dt = C* 

a n * 

(6) 


Eq. 6 is the crack growth relation used in this analysis. Taking the 
logarithm of Eq. 6a straight line relation, whose slope is n* and intercept 
is log C*, is obtained. 


log (da/dt) = log C* + n* log a (7) 

Several methods were tried in attempts to evaluate the parameters of 
Eq. 6. To determine the crack growth constants C* and n* for each growing 
crack individually the average crack growth rate between two inspections 
is calculated as the difference in crack sizes divided by the inspection 
interval evaluated at the average of the two crack sizes and at the mid 
point of the interval. Because two constants are to be determined a 
minimum of three successive inspections during which a crack grew are 
required to produce two points for Eq. 7 (Figs. 2 & 3). Not many cracks 
were permitted to grow through three inspections. 
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Due to this limitation only a relatively small number of the volu- 
minous data could be utilized to obtain crack growth parameters for 
individual cracks. 

2. Linear Least Squares Regression 

In order to increase the number of available data points an attempt 
was made to evaluate crack growth rate versus crack length data for each 
pair of successive inspections and then to use all of these values at a 
specific location on the aircraft for individual commands in a linear 
regression analysis as shown in Fig. 4. 

In a first attempt all the da/dt vs. a data were utilized in 

r mean 

order to find characteristic n* and C* values for each location. A 
computer program (Statistical Analysis System was used to perform a 
linear least squares fit to Eq. 7. The slope and intercept of the line 
are the n* and C* values, respectively. 

It was found that the exponent, n*, was lower than what had been 
anticipated from the examination of aircraft aluminum crack growth curves 
[11] and varied considerably. An n* of approximately 1.5 (n=3) had been 
expected. It is recognized that n* should be nearly constant in value 
since similar material and geometries are used in each aircraft and loca- 
tion. Referring to Eq. 4, C* varies with material constant C, material 
constant ji, geometric correction factor f(l/r) and stress range Aa. 

The variations in n* may be explained as follows: Fig. 5 represents 

two ideally growing cracks in the same location but on two different air- 
craft within the same command. The only difference between the two is 
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that aircraft B has experienced higher stress ranges than aircraft A. 

Note that both lines have approximately the same slope, n*, but differ in 
C* (Table 1) reflecting the difference in stress range. On the other 
hand a linear least squares fit of all data points (six) from both air- 
craft results in an n* = .92 instead of an average value of 1.01. If 
each crack is treated individually and average values of n* and C* are 
calculated less error is introduced. 

3. Non-linear Regression 

In another attempt to use all available data directly the rate of 
growth equation (Eq. 6) was integrated 

da/dt = C*a n * (6) 

therefore a” n da = C* dt 
Integrating both sides 


a" n * da = C* 


dt 


*1 


|V (a' 1 ’"*’ - a,.' 1 '"*’) ■ C (t - 


* 1 > 


,( 1_n * )= C* ( 1 - n* ) (t - t.) + a.^ 1 "^ 


a - 


C*(l -n*) (t - t.) + a, (1 " n * ) 


i ' i 


1-n* 


(8) 


where t. = time to crack initiation 

a. = crack length at crack initiation = .03 in. is obtained (Fig. 6). 
By assuming certain values of n* Eq. 8 may be simplified: 
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for n* = 2 


?■ - C * + 57 


where C-j = -C* 


For n* = 1 . 5 


1 

a 


= C-jt + C 2 





) + 



carrying out the operations and factoring yields 


( 9 ) 


where 


a “ C-jt 2 + C 2 t + C 3 


C, = 


C* c . c 
4 '* 4 


C* C^_ . 

y/i T~ 2 1 


00) 


C 3 = 


C* . 2 

~r 


/aT a i 


Equations 9 and 10 were evaluated again using regression analysis for 
all a^ versus .t data for a given location in each command (Fig. 6). Results 
were again unsatisfactory because n* values were afflicted with a large 
scatter. 
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4. Conjugate Gradients Method 

The next approach was an attempt to minimize the square difference 
between observed and predicted crack length by the method of conjugate 
gradients. 

The function to be minimized may be written as 


F = z (A.-a. ) 2 
k=l K K 



+ C* (1-n*) 


( V t i ) 


1-n* 


( 11 ) 


where = observed crack length 

a^ = predicted crack length 

t^ = observed time 

To minimize the function, its partial derivatives 


3F _ 
3C*" 


N 

Z 

k=l 


-2 (t k 




02 ) 



N 

z 

k=l 



(13) 


9F _ ? 2n* 

3n* k f 1 (1-n*) 2 


W - 0 ^] 


(i~*)l (t^f) 


In B, 


W'"" 


'In a, 


C*(t k - t 1 )] (14) 

are set equal to zero and the nonlinear equations, Eqs. 12-14 are solved 
simultaneously. 
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In the above relations 


B k = a. 0 '"*) + C*(l-n*) (t k -t.) 


The gradients, ^*, and , were used in a computer program (SSP 
subroutine FMCG), It was not possible to set bounding values in the sub- 
routine. Therefore, certain physical limitations, such as positive crack 

<rk> 

length, were violated. The term, B k , represents crack length and was 
in some cases forced into the negative region by the FMGC program in order 
to minimize the function. The method was consequently abandoned. 


5. Variable n* Method 

Crack growth parameters for individual cracks were also calculated and 
were then averaged to find representative C* and n* values. 

For this method only those da/dt versus a data were used in which 

a crack grew through at least two inspections (three points on an a^ versus 
t^ curve). Since stress ranges could vary greatly from one inspection to 
another, not all n* and C* values calculated were acceptable. Therefore, 
limitations were set on n* (1.0 < n* < 2.5), and on C* (.05 > C*). Many 
of the cracks that grew sufficiently to be investigated yielded n* and C* 
outside the acceptable range. 

The acceptable parameters for cracks within a given location and 
command were averaged to yield mean values of n* and C*. These average 
parameters were then used to extrapolate all cracks that grew, from their 
first sightings backwards to the time of crack initiation (t-j) and to 
initial flaw size (au). 

Table 2 lists the raw data of crack length and inspection times for 
the relevant locations. 
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This method yielded mean values of n* from approximately 1.2 to 1.9 
(Table 3). When these n* values were used to calculate a Q , the initial 
flaw size^ varied with n*; a lower n* produced a smaller a^ (see Figs, 7 
and 8). Since initial flaw size in an aircraft should be independent of 
the life of the aircraft, the values of a^ at various locations should be 
of the same order of magnitude. These results indicated that n* should 
indeed be a fixed material property. 

6. Constant n* Method 

A value for the material property n was selected after evaluating 
crack growth data on 7075 Aluminum from the Rockwell International Frac- 
ture Mechanics Data Bank [9] and comparing this to a NASA Langley study 
on laboratory crack growth of C-130 components [10], Based on these 
studies, n=3, i.e., n*=1.5, was chosen. 

With n* established, only one additional parameter needs to be deter- 
mined, hence individual crack growth data, da/dt versus a mean could be 
used (Fig. 9). It has been mentioned previously that when n* and C* were 
both allowed to vary, the values of n* and C* were frequently outside an 
acceptable range. By fixing n* and letting only C* vary, the results yielded 
a reasonable C* value every time. Due to the increase in the number of 
reasonable C*’s for all crack locations, more crack locations could be 
evaluated. 

The values of C* for the various growing cracks within a location were 
averaged to yield a C*. This C* value and n*=1.5 were used to extrapolate 
the a^ versus It data of the first crack that grew at a particular location 
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on an aircraft. The resulting 
It might be noted that all the 
as expected, of the same order 


values of t^ and a Q are shown in Table 3 
a^ values for the various locations were 
of magnitude (see Section III. 2). 
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Section III 


STATISTICAL DESCRIPTION OF THE DATA 

Statistical distribution functions have been fitted to the times to 
crack initiation and to initial flaw sizes. The probablity of time to 
crack initiation will be discussed first. 


1 . Time to Crack Initiation Distributions 
a. Two Parameter Wei bull Distribution 

Time to crack initiation (or failure) has in the past been described 
in terms of the Weibull distribution [24]. The two parameter Wei bull proba- 
bility of failure is given as 


F. (x) = 1 - e 
■i 


(x/B ) 01 


1.2 (a/x) is the shape parameter 
x/ r (1 + 1/a) is the scale parameter 


x = mean 

a - standard deviation 

This distribution fits most of the data but fails to fit the end 
points (see Figs. 10 through 17 and Table 4). The lower portion of the 
data (earliest crack initiation) is considered extremely important, there- 
fore it is necessary to try other statistical models that might fit these 
points better. 

b. Three Parameter Weibull Distribution 
The three parameter Weibull distribution 
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(16) 


F. (x) = 1 - e 



a 


has also been used. Because of the possible existence of flaws at time 
t=0 this distribution function will have a negative lower limit, x Q and 
Eq. 16 will have to be restricted to values for t^ >0. Consequently the 
distribution of crack initation times becomes: 


and 


F t (x) 

F t (0) 
L i 


F t (x) 
t i 



for x < 0 


for x = 0 


for x > 0 


(17a) 


(17b) 


(17c) 


that is, the probability that the crack initiation time is zero has a 
fixed, nonzero value. The corresponding density function consists of a 
Dirac function at t^ = 0 with an area given by Eq. 17b and a truncated 
density function £21], These functions are schematically shown in Figs. 

18 and 19, 

Such distributions have also been fitted to the data. There is 
significant improvement over the fit of the two parameter function, but for 
very short initiation times this distribution does not exhibit a good fit 
either. The distributions are presented in Figs. 20-27 and Table 5. 


c. Johnson Distributions 

Johnson distributions were examined next, due to the fact that the 
Johnson Sg distribution has four adjustable parameters [24]. In 
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order to determine which one of three types of Johnson distributions would 
be applicable, values of 3-, and 3 2 were calculated for each data set. 


3 


1 



(18) 


4 

(m 2 ) 2 


(19) 


where m 2 , m 3 , and m^ are second, third and fourth moments of the data, 
respectively: 


1 (x,-x ) 2 (20) 

2 n i=l 1 

nu = E (x.-x) 3 (21) 

6 n i=l 1 

% = ~ z (x.-x) 4 (22) 

* n i=l 1 


From an examination of Table 4 and Fig. 28 [24] it is noted that the 
Johnson or S B distributions apply. 

The Sg distribution has two shape parameters, y and n, one location 
parameter (lower bound), e, and one scale parameter, A. The distribution 
is bounded by e and e+A. In the case of time to crack initiation, the 
range. A, is unknown but can be estimated by the following equation: 

(x 0.5~ s) ( V e) * (x 0.5' e) (x l-o.~ e) ' 2 ( V e> (x l-cT e) 

,2 


A = (X 


0.5 


:) 


( X 0.5” E ^ " ( X cf e ^ ^ X l-a“ G ^ 


(23) 
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where x Q 5 = median of the data 
x q = a percentile of data 
x-j = (1 -a) percentile of data 

Since primary interest lies in fitting the extreme data points, a = 1% was 
used in most cases. More extreme values were often necessary to assure 
that x is positive. In some cases the term (x Q 5 -er - U a -e) ^ x i- a “ e ^ 
becomes negative if extreme values are not chosen. In the case of loca- 
tion 74-90 TAC, an x^ value had to be extrapolated below the lowest data 
point to obtain a positive x. 

Once e and X are established, the two shape parameters may be calcu- 
lated as follows: 



where Z = the stanrdard normal variate for a 100th percentile 

a 

Z n = the standard normal variate for (l L a) 100th percentile 
I -a 

Once e. A, n, and y have been calculated for each set of data, the 
following density function can be numerically integrated to find the 
distribution function 



( 26 ) 
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The results of the Johnson Sg distribution for time to crack initiation 
are shown in Figs. TO to 17 and the parameters for each location are presented 
in Table 4. 

Because x a and x^_ a may be arbitrarily chosen this family of distribution 
functions can be made to fit the short initiation time data points, but at 
the expense of missing the large number of data points in the central region. 

d. Statistical Distribution Based on Service Load and Fatigue 
Performance 

Because none of the afore discussed distribution functions fit the 
data equally well in the central region and the extremes, a technique based 
essentially on the concept of stress-strength interference will also be 
examined [6]. 

It is assumed that the fatigue strength, R, of the structure in a 
laboratory test can be described by a three parameter Wei bull distribution. 

The probability that the fatigue life is shorter than a given value, x is 
written as 

/ x ~ x o\ 

\e-x n / 

F r (x) = 1 - e 0 (27) 

where R and hence x is measured in laboratory hours. 

The statistical density function of the number of gusts that produce 
significant fatigue damage in the aircraft in t flying hours has been 
described by Yang [6] as a continuous version of the Poisson process 

f N (t)( x i) = (fct) 1 e" Xt /r(x 1 +l) 0 < x 1 < • (28) 

where the random variable, N(t), is the number of gusts in t flight hrs. 
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. X is the average number of significant gusts per flight hr, xt 
is the average number in t flight hrs and r(\) is the gamma function. 

Because the variance of the Poisson distribution is also xt, 
the coefficient of variation which is a measure of dispersion becomes 

V N(t) * < 29 > 

The number of gusts N(t) experienced by an aircraft in t flight 
hrs can be related to the number of equivalent laboratory hrs, S(t) as 

S(t) = N(t) As (30) 

where 5 is a socalled severity index that in the present analysis is 
assumed to be unity. As a consequence the laboratory spectrum is 
equivalent to the average flight spectrum. 

The density function, of gust loading, Eq. 28 can be transformed 
with the aid of Eq. 30 into a density function in terms of equivalent 
laboratory flight hours 

f S( t )(x) = Xt xX e' Xt r '(Yx+T)" (31) 

The probability that a crack will initiate in t flight hours is 
equivalent to the probability that the fatigue strength, R, is 
exceeded by the service loading converted to laboratory flight hrs, 

set) 
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F T (t) = P[R < S(t)] (32) 

Eq. 32 is a form of the classical stress-strength interference rela- 
tion [25] and can be expressed as 

F ? (t) = ( F r (x) f s(t) (x) dx (33) 

.0 

Where F R (x) is given by Eq. 27 and f s ( t ) ( x ) Ec, ‘ therefore 



For convenience of integration a change of variables is introduced: 

y = ax , dy = xdx (35) 

Substituting into Eq. 34 the probability of failure becomes 


a 



Eq. 36 is a convolution integral that requires numerical integration. 
The values of ct, $ and x Q are estimated by fitting a three parameter 
Wei bull distribution to the central part of the crack initiation time 
data (Table 5). 

Because l//xg" is a measure of the dispersion in g hrs, a small 
value of a produces a large dispersion. Choosing x values by trial 
and error (xg is usually less than 100) Eq.36 is evaluated and is 
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fitted to the crack initiation time data so that the curve approaches 
the data in the lower tail (short initiation times) of the distribution 
as shown in Figs. 29 to 36, 

Smaller values of x tend to bend the lower tail upward and the 
upper tail downward. 

Because short crack initiation times are of interest, the poorer 
fit at the upper tail is not considered significant. 

2. Initial Flaw Size Distributions 

It is not unexpected that material imperfections, manufacturing 
methods and the limitations of inspection techniques produce a certain 
number of undetected minute cracks even before the aircraft is put into 
service. Once the average parameters of the crack growth relation are 
evaluated and t^ is calculated for each crack, Eq. 8 can be further 
extrapolated backwards to time t=0 to yield an initial crack size as 
shown in Fig. 2. Only data with extrapolated initial crack sizes 
greater than zero can be used. 

Initial flaw sizes were plotted on normal probability paper. The 
plots indicated that the distribution functions were not normal; the 
data did not fall on straight lines. As an alternative the Johnson 
Sg distribution was again employed. Since .03 inches (.0762 cm) 
had been chosen as the crack initiation size, the nondefect population 
(those locations with initial flaws less than .03 inches) ranges from 0.0 
to 0,03 inches. Thus, e = 0.0 and X = 0.03. The shape parameters n and 
x, were found as previously shown from Eqs. 24 and 25 with a = 1-2%. 
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The parameters are listed in Table 6 . Distribution functions, are 
presented in Figs. 37 to 44. it is seen that the Johnson S g distri- 
bution fits these extrapolated points very well. Because initial 
flaw sizes should be independent of loading their distributions 
should be similar. It is seen from Table 6 that the means and 
standard deviations are of similar magnitude. As a consequence all 
initial crack sizes could be considered to belong to the same 
population. An overall distribution function has been derived for 
the combined data and is presented in Fig. 45. The mean initial 
crack size is .00432 in. and has a coefficient of variation of .2%. 

The two Commands considered here, TAC (Tactical Air Command and 
PACAF (Pacific Air Command), had 263 and 223 C-130 Aircraft inspected 
respectively. 

Since two symmetric locations and two holes at each location 
were all considered to be identical the total number of potential crack 
sites at one location is 1052 and 892 for TAC and PACAF respectively. 

Of these 16% in TAC and 30% in PACAF have actually developed growing 
cracks at the most vulnerable location (75-91, see Table 3)- Aircraft 
operated by PACAF have a greater proportion of growing cracks at all 
locations while aircraft in the Tactical Air Command developed fewer 
inspectable and hence extrapolatable cracks. This fact would indicate 
that PACAF aircraft have been subjected to more severe loading condi- 
tions. Such a conclusion is not always born out by the data shown in 
Table 3. The value of C* for locations 75-91 is greater for PACAF 
than for TAC but at the same time other locations indicate an opposite 
trend. 
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Section IV 


THE RELIABILITY OF INSPECTION 

1 . Distributions of Detected and Undetected Cracks 

The probability that an existing crack is detected is a function of a 
number of circumstances. These include: a.) The method of detection such 

as radiography, thermography, dye penetrant, ultrasonic or visual inspection* 
etc., and b. ) The sensitivity and accuracy of the instrument or the alert- 
ness of the inspector. These conditions are influenced primarily by the 
size of the crack at the time of inspection. 

As a consequence not all cracks are found and the probability, P[D|a], 
that a crack is detected, given that its length is a^, is a function of the 
crack size. 

The methods of inspection or the reliability of the procedure used for 
crack detection in 0130 airplanes are not available. Data can however be 
found in the literature on the probability of crack detection in aluminum 
using various ultrasonic techniques [26]. A set of data from Ref. 26 is 
shown in Fig. 46 with a Weibull type probability function 



where "a = 2.0 and "e = .075 in. (.19 cm), fitted to it. Similar proba- 
bility functions have been incorporated in wing failure analyses and have 
been reviewed by Eggwerts [27]. 

The C-130 inspection data indicates that cracks shorter than .03 in 
have not been found. The assumed probability function used here predicts 
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that a small fraction, 15%, of such cracks are located, while the probability 
of detecting a crack .2 in. in length is almost unity. 

If the distribution of crack sizes prior to inspection is known, the 
distributions of detected and undetected cracks may be calculated utilizing 
the probability of detection function, Eq. 37 in conjunction with Bayes 
theorem [28,29], 

The probability that a flaw of size, a^, will be between x and x+dx 
given that the flaw is detected may be written as 


P[x<a<(x+dx) | D] = 


P[D|x<a<(x+dx)]P[x<a<[x+dx] 


where D stands for detection, P[D|x<a<(x+dx)] is the probability of detecting 
a flaw of size a_ (from Eq. 37), P[x<a<(x+dx)] is the probability that the 
crack size is a^ and P[D] is the total probability of detection independently 
of size. Eq. 38 may be expressed in terms of density functions as 

f a [x|D] = k ] P[D | a] f a (x) (39) 

in which f, (x) is the density function of flaw sizes prior to detection, 

a 

f a [x|D] is the density function of detected flaw sizes (posterior density), 
P[D | a] defined earlier, is called the likelihood function in Bayes' theorem 
terminology and ky= 1/P[D] is a normalizing constant that will make the 
area under the f a (x|D) curve equal to unity. 

In a similar manner the density function of undetected cracks may 
also be derived: 

f a [x|D] = k 2 P[D|x] f (x) (40) 
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where D stands for no detection, P[cT| x] = 1 - P[D| x], and k 2 is a second 
normalizing constant. The above technique has been utilized by Davidson 
[21] in a slightly different form with exponential detectability and 
crack size distribution functions assumed. 

Adding Eqs. 39 and 40 


f a (x|D) 

— + 

■1 



C41) 


indicates that the reciprocal normalizing constants, l/k-j and l/k 2 represent 
the fractions of detected and undetected cracks. 

Because the density function of cracks before inspection is usually 
unknown, the posterior function, f a £x | D] , of detected cracks may be used 
to evaluate the prior, f (x),and the density of undetected cracks* 

3 

f a tx | , from Eqs. 39 and 40 [29]. 

As indicated by Eq. 8 crack size is a function of the parameters a. 
and t. and the time of inspection t. As a result Eqs. 39-40 are also 
dependent on these quantities. Since inspections were carried out at ir- 
regular intervals the above relations could only be utilized if detected 
cracks were mathematically grown to a common inspection time. Alterna- 
tively Eqs. 38-41 can be transformed to yield density functions of crack 
initiation times, t.., at a common crack size, a_. = .03 in (.76 cm). 


2. Distributions of Detected and Undetected Crack Initiation Times 

Eqs. 39 and 40 may be rewritten in terms of crack initiation times, 


f [x|D] - k 3 PtDjt,] 


f t. (X> 


(42) 
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and 


f t .[x|D] = k 4 (l-PCD|t t 3)f t (x) 


(43) 


where f^[x|D] and f^.[x|D] are the density functions of detected and un- 
detected crack initiation times, P[D 1 1^ 3 is the likelihood function for 
detection and f t . (x) is the prior density function of all crack initiation 
times. 

P[D|t..] is derived from Eq. 37 by substitution of Eq. 8. Hence 


P[D|t..] = 1 - exp 



f C*(l-n*)(t-t 1 ) 




It should be noted that Eq. 44 becomes equal to unity when 


(44) 


, (1-n*) 
(t " V = ' C*(l -n*) 


(45) 


because within this time interval a crack will grow to infinite length. 
Hence P[D|t..] will be equal to Eq. 44 for ( t-t . ) < -[a.. ^~ n * Vc* 0 -n*)] 
and becomes equal to unity for ( t-t .) >_ -[a.. ^~ n *Vc*(l-n*)], which 
indicates that the probability of detecting an initiation time, t.., 
depends on the time of inspection, t. As a consequence the density func- 
tions of detected and undetected initiation times will also become func- 
tions of the inspection time while the density function of all initiation 
times should be independent of when the inspection is carried out. 

The distributions of t^ described in Section II were derived from 
long-term records, in some instances covering 7000 hrs of flying time. 
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It is therefore not unreasonable to assume that the great majority of cracks 
have grown to inspectable size and that distributions of crack initiation 
times derived from these by backward extrapolation are good approximations 
of the total initiation time distributions. Such an assumption is conserva- 
tive because long initiation times (small cracks) make up the majority of 
still undetected initiations. 

With this assumption it is possible to determine the density functions 
of discovered and undiscovered initiation times for an inspection carried 
out earlier in the life of the structure. 

The three parameter Wei bull distributions, Eqs. 17, fitted to the 
derived data of Section II (the convolution integral representation is 
mathematically more complicated) 



will be used where <5(x), is a Dirac function. 

For specified inspection time, t, crack growth parameters C* and h* 
and a. = .03 in., Eqs. 44 and 46 substituted into Eqs. 42 and 43 yield the 
density functions of detected and undetected initiation times. As in the 
case of detected and undetected cracks the reciprocal normalizing coeffi- 
cients add to unity 



(47) 


and indicate the fractions of discovered and undiscovered initiation 
times. 
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The results of such calculations are shown at location 76-92, TAC 
Command in Figs, 47-49 for inspection times t = 2500,4000 , and 5000 
hrs. The parameters of the distribution functions are presented in Table 7. 


It should be noted that the density functions for detected and un- 
detected initiation times have been divided by their respective normalizing 
constants so that the areas under the two curves are proportional to the 
fractions of discovered and undiscovered initiation times and the two areas 
sum up to the area of the density function for all initiation times. 

An examination of Figs. 47-49 indicates an expected trend. If 
inspection is performed too early, the majority of initiation times will 
remain undetected while more and more crack initiation times are found at 
later inspections. 


The aforegoing discussion has utilized longterm inspection data to 

derive the prior distribution of crack initiation times. Such records are 

however usually unavailable for aircraft in service. On the contrary at a 

particular inspection time the distribution of observed crack initiation 

times is actually the posterior distribution, f t (x|D), while the prior 

■ i 

density, f. (x), of all initiation times and the density function of undetected 
initiation times, f. (x|D), are unknown. 

Because the latter two functions are of primary interest Eqs. 42 and 43 


may be rewritten using Eq. 47 as 

f t [x|D] 

f t (x) ■ __i (48) 

1 k 3 P[D| t 1 ] 


and 


f t (x |'D) 
u i 


(l-P[D|t.])f. [x|D] 
1 r i 

(k 3 -l) P[D|t.] 


(49) 
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These equations would predict the density functions of all crack initiation 
times as well as that of the undetected times provided that inspections 
are carried out at approximately the same number of flight hours on all 
aircraft in a fleet. 
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Section V 
CONCLUSIONS 

The methods presented indicate the feasibility of determining the 
statistical distributions of initial crack sizes and times to crack 
initiation using backward extrapolation of inspection data. It appears 
that of the two crack propagation parameters the value of n* is a 
material constant with a value of n* = 1.5 while the second parameter C*, 
in addition to being dependent on material and geometry, is also a 
function of the stress range. 

Evaluation of C* for each individual growing crack and averaging 
such values for specific crack locations and military commands yields 
satisfactory distributions of crack initiation times and initial crack 
sizes. 

Crack initiation times are dependent on both the initial crack 
size and on C*, the latter being dependent on location and load, while 
the distribution of initial crack sizes is independent of these factors. 

The three parameter Wei bull distribution fits the crack initiation 
times reasonably well. The convolution integral, that involves the 
three parameter Weibull distribution of fatigue strength and the Poisson 
distribution of gust occurrences does not improve the fit appreciably 
because this method has been developed for the purpose of evaluating 
laboratory fatigue tests under simulated load spectra rather than for 
the analysis of service data. Two parameter Weibull or Johnson S B 
distributions do not fit the data very well. 
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The distributions of initial crack sizes on the other hand are 
well represented by the Johnson Sg distribution. The methods presented 
permit evaluation of a single initial crack size distribution for all 
growing cracks independently of the different loads in the two commands 
or the location of cracks on the wing. Though only four locations were 
analyzed it is believed that the distribution of initial crack sizes is 
representative of other locations as well. The information can con- 
sequently be used for the prediction of crack growth for the C-130 air- 
craft as well as for structures made of the same material with similar 
geometry and having undergone similar manufacturing processes. It is 
recognized that the distributions of time to crack initiation and initial 
crack size are not independent of each other but are interrelated through 
Eq. 8. Consequently a probability transformation can be performed on 
one of the two distributions to derive the density function of the other. 
It is obvious that such a transformation will result in distribution 
functions different from those presented here though it is expected that 
the derived function would be well approximated by the fitted curves. 

A transformation will be carried out at a future date. 

The Bayesian approach presented here illustrates a method for esti- 
mating the proportion of undiscovered crack initiation times at various 
inspections and indicates the influence of the reliability of crack 
detection methodology particularly during early inspections when cracks 
are small . 

Because the information concerning the reliability of the particular 
method of nondestructive examination used for the C-130 aircraft was not 
available an assumed likelihood function was used. For future systems 
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the documentation of NDE reliability should become an integral part of 
service inspection. 

A similar approach may be used to determine the fraction of unde- 
tected crack sizes at a particular inspection time. For this purpose 
cracks would have to be analytically grown to common inspection times 
and their distributions analyzed. Such work may also be carried out 
in the future. 

The present analyses have utilized average values of the crack 
propagation parameter C*. The statistical distributions of the param- 
eter have not been evaluated. Again calculations, with C* as a random 
variable, may be performed at a later date. 
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Crack Growth Constants for Two Hypothetical Aircraft 



36 



TABLE 2 


C-130 Inspection Data for Selected Locations 
(Continued p. 38-60) 


Explanation of columns 1-8 


Column 1 . 
Column 2. 
Column 3. 
Column 4. 
Column 5. 
Column 6. 

Column 7. 
Column 8. 


Aircraft Identification Number 
Location of Rivet Hole (see Fig. 1) 

Two Holes at Each Location: (I) Inboard and (0) Outboard 

Direction of Crack Growth: (F) Forward and (A) Aft 

Crack No. at Rivet Hole. Only First Growing Crack Used 

User Command (P) Pacific Air Command, PACAF; (T) Tactical 
Air Command, TAC 

Length of Crack in Hundredths of Inches 
Inspection Time in Flight Hours 
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Crack Propagation Parameters for Constant 
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al Parameters for Time to Crack Initiation 
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Statistical Parameters for Fitting Time to Crack Initiation 
Three Parameter Wei bull Distribution and Convolution 
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Statistical Parameters for Initial Flaw Sizes 
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TABLE 7 


Normalizing Constants and Percentage of Detected and 
Undetected Crack Initiation Times 


Inspection Time 
t hours 

k 3 

% 

Detected 

k 4 

% 

Undetected 

2500 

16.60 

6.02 

1.06 

93.98 

4000 

1.60 

62.56 

2.67 

37.44 

5000 

1.01 

98.56 

69.34 

1.44 
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76-92 73-89 



75-91 *- 74-90 




CRACK LENGTH 



INSPECTION TIME, t, (hrs) 


Fig. 2. Extrapolation of Inspection Data. 
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Fig. 3. Determination of Crack Growth Parameters 
Variable n* Method. 
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Fig. 4. Determination of Crack Growth Parameters. 
Linear Regression Method. 
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CRACK GROWTH RATE (da/dt) 
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Fig. 6. Determination of Crack Growth Parameters. 
Nonlinear Regression Method. 
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PROBABILITY DENSITY - f(a n ) 



inches 


0 .01 .02 .03 .04 .05 

centimeters 

INITIAL FLAW SIZE - a Q 

Fig. 7. Histogram of Initial Flaw Size Using a 
Larger n*. 
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Fig. 9. Determination of Crack Growth Parameters. 
Constant n* Method. 
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PROBABILITY 



Fig. 10. Distribution of Crack Initiation Times • Two Parameter 
Wei bull and Johnson Sg Distributions. 
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PROBABILITY 



Fig, n. Distribution of Crack Initiation Times. Two Parameter 
Weibull and Johnson Sg Distributions. 
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PROBABILITY 



Fig. 12. Distribution of Crack Initiation Times. Two Parameter 
Wei bull and Johnson S B Distributions. 
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PROBABILITY 



Fig, 13. Distribution of Crack Initiation Times. Two Parameter 
Weibull and Johnson S B Distributions. 
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PROBABILITY 



Fig. 14. Distribution of Crack Initiation Times. Two Parameter 
Weibull and Johnson Sg Distributions. 
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PROBABILITY 



TIME TO CRACK INITIATION, t., (TO 3 flight hours) 

Fig. 15. Distribution of Crack Initiation Times. Two Parameter 
Wei bull and Johnson S B Distributions. 
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PROBABILITY 



Fig. 16. Distribution of Crack Initiation Times. Two Parameter 
Wei bull and Johnson Sg Distributions. 
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Fig. 17 . Distribution of Time to Crack Initiation. Two Parameter 
Wei bull and Johnson Sg Distributions. 
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0 2 4 6 8 10 


TIME TO CRACK INITIATION, t., (TO 3 flight hours) 
Fig. 20. Distribution of Crack Initiation Times. 
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PROBABILITY, Fj.(t.) 



0 2 4 6 8 10 


TIME TO CRACK INITIATION, t i , (TO 3 flight hours) 
Fig. 22. Distribution of Crack Initiation Times. 
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PROBABILITY 



0 2 4 6 8 10 


TIME TO CRACK INITIATION, tj, (10^ flight hours) 
Fig. 23.. Distribution of Crack Initiation Times. 
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PROBABILITY, F x .(t.) 



TIME TO CRACK INITIATION, t i , (TO 3 flight hours) 
Fig. 25. Distribution of Crack Initiation Times. 
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PROBABILITY 



TIME TO CRACK INITIATION, t 4 , (TO 3 flight hours) 
Fig. 26. Distribution of Crack Initiation Times. 
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PROBABILITY 



TIME TO CRACK INITIATION, t i , ( 1 O 3 flight hours) 
Fig. 27. Distribution of Crack Initiation Times. 
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Fig. 28. Johnson Distribution 
Applicability Chart. 
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PROBABILITY 



TIME TO CRACK INITIATION, t. , (10 3 flight hours) 
Fig. 29. Lower Tail of Crack Initiation Time Distribution 
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TIME TO CRACK INITIATION, tj, (TO 3 flight hours) 

Fig. 30. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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Data Points 



0 12 3 4 

TIME TO CRACK INITIATION, t 1 , (10 3 flight hours) 


Fig. 32. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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PROBABILITY, 


74-90 TAC 


— Convolution 
• Data Points 



0 1 2 3 4 


TIME TO CRACK INITIATION., t 1§ (10 3 flight hours) 

Fig. 33. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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Fig. 34. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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0 12 3 ^ 

TIME TO CRACK INITIATION, t is (TO 3 flight hours) 


Fig. 35. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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PROBABILITY, 



0 12 3 4 

TIME TO CRACK INITIATION, t-j , (10 3 flight hours) 

Fig, 36. Lower Tail of Crack Initiation Time Dis- 
tribution. 
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PROBABILITY 



Fig. 37. Distribution of Initial Flaw Sizes* 
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PROBABILITY 



Fig. 38. Distribution of Initial Flaw Sizes • 
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PROBABILITY 



Fig. 39. Distribution of Initial Flaw Sizes. 
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PROBABILITY 



Fig. 40. Distribution of Initial Flaw Sizes. 
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PROBABILITY 



Fig. 41. Distribution of Initial Flaw Sizes. 
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PROBABILITY 



Fig. 42 . Distribution of Initial Flaw Sizes. 


107 



PROBABILITY 



Fig. 43. Distribution of Initial Flaw Sizes. 
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Fig. 45. Distribution of Initial Flaw Sizes for all 1256 
Growing Cracks • 
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Fig. 47. Distributions of Crack Initiation Times for a 
Typical Wing Location for an Early Inspection. 
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PROBABILITY DENSITY - f (t- ) 



CRACK INITIATION TIME - 

Fig. 49. Distributions of Crack Initiation Times for a 


Typical Wing Location for a Late Inspection. 
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